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TECHNICAL PAPER 


HYPERSONIC RAREFIED WAKE CHARACTERIZATION 

INTRODUCTION 


Due to the aeroassist flight experiment (AFE), there has been a renewed interest in rarefied 
hypersonic flight after a considerable period of dormancy. The AFE program, which is now on hold, was 
a mission to obtain science data at rarefied flow conditions not available to the aerothermodynamicist in 
wind tunnels. Those science data are needed to confidently design an aeroassisted space transfer vehicle 
(ASTV). 

An ASTV will consist of an aerobrake, a payload, and a propulsion unit. The aerobrake protects 
the payload from the severe aerothermodynamic environment of hypersonic entry into the atmosphere. 
However, the environment of the payload which resides in the wake of the aerobrake still must be 
defined. Little is known about these wakes. It is the purpose of this paper to characterize the rarefied 
hypersonic wake and its interaction with a payload. This basic information will be useful in designing 
the aerobrake and its payload and for understanding the tethered satellite wake. 


COMPUTATIONAL METHOD 


The direct simulation Monte Carlo (DSMC) method 1 ” 3 is the only method presently available to 
handle rarefied wake flow. The five-specie air chemistry model 4 was used in these calculations. The 11- 
specie model, which includes the five species O 2 , O, N 2 , N, NO and their ions along with the electron, is 
appropriate for these reentry conditions, but lacks the maturity of the five-specie model. It is hoped that 
the ionization which will occur will have little influence on the structure of the wake. Dr. Bird’s two- 
dimensional or axisymmetric code 5 was used for the axisymmetric calculations. An equivalent code 
written by the author was used for the three-dimensional AFE calculations. 


RESULTS AND DISCUSSION 


While the capabilities and accuracy of the DSMC method are quite well documented, benchmark 
calculations are still necessary to add to the experience base when it is used to predict results for new 
problems or different geometries. While Dogra et al. 7 looked at hypersonic rarefied flow past spheres, 
including the wake structure, the results were compared only with the drag force from Legge and 
Koppenwallner’s experiment. 8 


Sphere Wake Flow Benchmark 

Muthoo and Brundin 6 presented, at the RGD symposium in Gottingen, near wake measurements 
of a sphere at M = 5.6 and Kn = 0.005 run in the Oxford University low-density wind tunnel. These data 


were used as a benchmark case in order to further validate the code for this study. A cold gas model used 
T’ref = 300 K and viscosity exponent co = 0.91 rather than oj = 0.75 which is normally used for air at 
normal temperatures because of the very cold freestream temperature of the Oxford tunnel. 

The DSMC results are given in figures 1 through 3. Contour plots of density, velocity, and tem- 
perature are shown in figure 1. The wake density field from both analysis and experiment indicated a 
broad recompression region and weak trailing shock or recompression fan. There is no recirculation in 
the base region as indicated by the vector plot given in figure 2. Experiment and DSMC are in excellent 
agreement with respect to the defects of velocity and density along the wake centerline beyond an x/d of 
one as shown in figure 3. The translational kinetic temperature was calculated from the collection of 
molecules in each cell, and the pressure was calculated from the equation of state in order to compare 
with the experimental data. For regions of nonequilibrium, the thermodynamic temperature and the 
translational temperature will not be equal. Both analysis and experiment (fig. 3) show that there is no 
separation in the wake; however, in their paper, Muthoo and Bruden state that rear stagnation point is 
located approximately at x/d - 1.5 where the axial static pressure is a maximum. The DSMC results 
indicate that the maximum static pressure occurs near an x/d =1.25; however, an interpretation that this 
is a stagnation point is not correct. This is a characteristic of rarefied wakes, and more will be said later. 
Obviously, experiment and analysis agree, as is shown in figure 3, as to the velocity distribution in the 
wake, and there is no flow reversal. 

The differences between the DSMC results and experiment in the near wake might be attributed 
to experimental error due to the large mean free path in that region of the flow. The physical size of the 
instrumentation, the low-density corrections required on their response, and probe interference with the 
flow probably significantly affect the experimental results in close proximity to the sphere. The flow in 
this region is not in equilibrium and with the rarefaction the continuum concepts of temperature and 
pressure become invalid. Density remains the significant flow property, but the experimental value was 
not measured directly, of course. It is a result of analysis of hot wire and pressure data. Even today, it 
would be very difficult to make meaningful measurements in the very near wake of a sphere. 


Sphere Wake Flow Characterization 

The surface shear stress near the sphere aft stagnation point indicates that the flow may be near- 
ing separation for the benchmark flow conditions. Using the same grid, the calculations were pushed to 
higher Reynolds numbers until flow separation occurred. Further increases in Reynolds number then 
increased the extent of the separation zone. Several cases were run at various Reynolds numbers and 
Mach numbers as indicated in tables 1, 2, and 3. Typical TJT 0 = 1 results are plotted in figures 4 and 5 
for Re 0 greater than 600. The calculations at the highest Reynolds numbers are certainly coarse grid 
results near the windward stagnation point, but on the wake side the grid is adequate relative to the local 
mean free path. However, near the sphere, the steep gradients indicate that a finer grid is desirable. In 
any case, the extent of the predicted flow separation is given in figure 6 as a function of the Reynolds 
number with viscosity evaluated at the flow total temperature, Re 0 . The value of Re 0 = 600 corresponds 
approximately to incipient separation for a TJT 0 = 1 for M = 5.6 and M = 30 as well, exhibiting some 
degree of Mach number independence. 

For a cold wall, one would expect separation to occur at lower Reynolds numbers, and that is 
what happens as shown in figure 6. Typical TJT 0 = 0.1 results are plotted in figures 7 and 8. While the 
data shown for the M = 30 case are insufficient to determine the separation Reynolds number, it appears 
that Mach number independence does extend to the cold wall case. 
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ASTV Configuration 


Since there are expected to be configuration effects upon the wake flow, a geometry more repre- 
sentative of an ASTV than a sphere has been considered. An ASTV will consist of an aerobrake and a 
payload. The configuration used here is a hemispherical aerobrake with a cylindrical payload for sim- 
plicity. The radius of the payload was arbitrarily set at half the radius of the hemisphere. The front sur- 
face of the aerobrake was set to 1,000 K, while the aft aerobrake surface and the payload were set to 
300 K. 


The freestream velocity was set at 10 km/s. The ambient conditions used correspond approxi- 
mately to the altitudes of 86, 92, 100, and 110 km and were taken from the 1972 Standard Atmosphere. 
Since the hemisphere diameter is 1.0 m, the Knudsen number is equal to the mean free path. The free- 
stream conditions are given in tables 4 and 5. 

For all the cases run over the transitional flow regime, the near wake was separated near the 
hemisphere. This is a marked contrast to the sphere results where the wake is not separated in the transi- 
tional flow regime. In fact, when molecular collisions were suppressed, i.e., collisionless flow, the wake 
still exhibited a separated flow streamline pattern. A representative wake flow field is shown in figures 9 
and 10 where density, velocity, and temperature contours are plotted. 

One might expect the flow over a hemisphere to separate at a lower Reynolds number relative to 
a sphere, but its persistence to free molecular flow indicates that the payload or afterbody must also be 
having an effect. One could eliminate the payload from the geometry to determine its effect; however, its 
effect can be eliminated by applying a plane of symmetry boundary condition at its surface. There is 
then no heat transfer or skin friction at the payload surface either. Thus, applying a specular reflecting 
boundary condition at the payload surface gives the advantage of keeping the same computational space 
and grid as before while eliminating the effects of the payload. Further, starting with the diffuse 
boundary condition results, the new results are obtained quickly. 

The results are given in figures 1 1 and 12. The wake is now unseparated (except at the lowest 
Knudsen number) in conformity with the sphere results. A comparison of the streamlines for the two 
cases is shown in figure 13. The location of the stagnation point on the payload is quite close to the 
aerobrake. In contrast to high Reynolds number flows, where this point is usually associated with high 
heating, the thermal environment at this point in the rarefied wake is quite benign. 

The heating, pressure, and skin friction distributions along the payload are shown in figures 14, 
15, and 16, respectively, for selected Knudsen numbers. The surface pressure peaks at a slightly lower 
x/d (0.6 to 0.8) value in comparison to the heating rate and shear stress (0.8 to 1.0). The heating at the 
stagnation point on the payload ( x/d < 0.15) is an order or so less than the maximum. 

AFE Configuration 

An ASTV will be a lifting body and will have a three-dimensional wake like the AFE configura- 
tion. The AFE payload is the carrier vehicle. The DSMC method was used previously to determine the 
wake flow characteristics of the AFE. 9 The wake results were a little puzzling, and one always questions 
the validity of first results. The AFE and its wake are three-dimensional and require a great deal of com- 
puter time and memory. The grid resolution is never as good as one would like with the computer 
resources being the constraint. In this case, the computations were done on the best facilities available, 
being the NAS facility at Ames Research Center. 
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Shown in figure 17 are the flow vectors obtained in the plane of symmetry of the AFE wake at an 
altitude of 100 km and a velocity of about 9.9 km/s. The flow separates from the aerobrake and impinges 
on the carrier vehicle (the afterbody or payload) near the upper part of the figure. Recirculation appears 
near the carrier vehicle in the lower part of the figure. The wake behind the afterbody is not separated. 
These wake characteristics were the same at all altitudes calculated from 95 to 1 20 km (Knudsen range 
from 0.03 to 0.8). This wake behavior is consistent with and now explainable by the ASTV results pre- 
sented in the previous section of this report. 


CONCLUSIONS 


_ The wake flow field has been calculated using DSMC for hot and cold spheres at freestream 
conditions near continuum. Separation on nearly adiabatic spheres has been predicted to occur near a 
postshock Reynolds number of 600 and for cold spheres near 100. For more rarefied conditions, the 
wake does not separate. It has been shown that an afterbody or payload in the wake of a hemispherical 
aerobrake will cause the near wake flow to separate over the whole transitional flow range from con- 
tinuum to free molecular flow. This is in sharp contrast to aerobrakes without afterbodies where the flow 
is generally not separated except in the near continuum. For a short, blunt-based afterbody like the 
carrier vehicle in the wake of the AFE aerobrake, separation occurs over the afterbody, while the wake 
behind the afterbody (carrier vehicle) is not separated. 

No matter how short the afterbody, the flow will be separated near the aerobrake end of the pay- 
load. However, the heating, shear stress, and pressure distributions will peak a considerable distance 
downstream of the attachment point. For example, for the ASTV case 2, the reattachment length is 
0.15 while the pressure peaks at an x/d of about 0.8 and the heating rate at an x/d of 1.0. It is clear that 
the heating on the payload can be sharply reduced by keeping the payload short. 

As previously mentioned with regard to the sphere benchmark case, the peak pressure in the 
rarefied wake does not correspond closely with a stagnation point as is the case with wakes at a high 
Reynolds number. The local flow is quite rarefied in the recirculation region with the mean-free-path 
being approximately equal to the reattachment length. In fact, the separation or recirculation which 
appears with the long time averaging of the molecular motion, where there are too few collisions for 
equilibrium to exist is really not a “separation” in the continuum sense, and the interpretation that this is 
a recirculation region stems from the continuum streamline concept. While this makes sense with respect 
to continuity (mass flow), it can give an erroneous impression of the actual local flow conditions. For 
example, consider a two-dimensional trailing edge problem with upper and lower streams expanding 
around the trailing comers and intersecting each other. In continuum flow, the two streams interact and 
turn to accommodate each other. In a very rarefied flow, the two streams may pass through each other 
with little or no interaction. The resulting streamline patterns from the two very different flows would 
look similar. The point being that, in addition to temperature and pressure, streamlines and consequently 
separation and reattachment do not have the same implication in rarefied flow that one is accustomed to 
in continuum flow. 
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Table 1. Freestream conditions for sphere calculations. 


Case 

V (m/s) 

n 

T (K) 

1 

721 

0.4068xl0 22 

41.25 

2 

721 

0.8137X10 22 

41.25 

3 

721 

0.2369x1 0 23 

41.25 

4 

721 

0.3353x1 0 23 

41.25 

5 

721 

0.5345x1 0 23 

41.25 

6 

2,827 

0.2935x1 0 23 

22.10 

7 

2,827 

0.6662x1 0 23 

22.10 


Table 2. Selected sphere freestream parameters. 


Case 

M 

M 

Re 

Re 0 

Kn 

5 

1 

28.97 

5.6 

627.0 

97.3 

0.0134 

4.69 

2 

28.97 

5.6 

1,254.0 

194.5 

0.0068 

4.69 

3 

28.97 

5.6 

3,652.0 

566.3 

0.0023 

4.69 

4 

28.97 

5.6 

5,169.0 

801.6 

0.0017 

4.69 

5 

28.97 

5.6 

8,240.0 

1,278.0 

0.0010 

4.69 

6 

28.97 

30.0 

33,058.0 

565.0 

0.0014 

25.1 

7 

28.97 

30.0 

74,600.0 

1,272.0 

0.0006 

25.1 


Table 3. Molecular parameters. 


Case 

CO 

TVefCK) 


Do 2 (mx 10 10 ) 

1-5 


300.0 

4.374 

4.203 

6-7 


1,000.0 

3.405 

3.265 


Table 4. Freestream conditions for ASTV calculations. 


Case 

V (km/s) 

n 

T (K) 

Mass Fractions 

Xch Xn 2 Xo 

1 

10 

0.1350xl0 21 

186.87 

m 


0.0000 

2 

10 

0.4640xl0 20 

186.87 

■ mm 

1 

0.0000 

3 

10 

0.1189xl0 2 ° 

195.00 

mm 


0.0365 

4 

10 

0.2144xl0 19 

240.00 

mmm 

HB9 

0.1080 
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Table 5. Selected ASTV freestream parameters. 
























(a) Nondimensional density, p/p^ (b ) Nondimensional velocity, V/V 



x/d 

(c) Nondimensional Temp, T/T D 

Figure 4. Selected variable contours about a sphere, case 4, TJT 0 = 1. 



0.5 0.6 0.7 0.8 

x/d 


Figure 5. Flow vectors in a near wake of a sphere, case 4, TJT 0 - 1. 



■ , M=30; Closed sym T w /T o =1.0; Open sym 
T w rr o =0.i 


Figure 6. Extent of laminar flow separation behind a sphere from DSMC. 




z/d 

(a) Nondiroensional density, p/ 



x/d 

(b) Nondimensional Velocity, V/V „ 



x/d 


(c) Nondimensinal Temp, T/T 0 

Figure 7. Selected variable contours about a sphere, case 4, T w /T 0 =0.1. 
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Figure 10. Flow vectors in an ASTV wake with diffuse B.C. on payload, case 2. 



(a) Nondimensional density, p/p 



(b) Nondimensional velocity, V/V w 



Figure 1 1. Selected contours in an ASTV wake with specular B.C. on payload, case 2. 
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(a) Case 1 



x/d 


(b) Case 2 




(c) Case 3 (d) Case 4 


Figure 14. Nondimensional heating distribution along payload. 
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(d) Case 4 


Figure 15. Pressure coefficient distribution along payload. 
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